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A  percolation  theory  based  model  considering  particle  size  and  its  distribution  is  proposed  for  composite 
electrodes  of  solid  oxide  fuel  cells  (SOFCs).  The  model  calculation  agrees  excellently  with  3D  numer¬ 
ical  reconstruction  results,  suggesting  great  validity  of  prediction.  Moreover,  it  is  also  consistent  well 
with  experiment  for  real  LSM  (lanthanum  strontium  manganite)-YSZ  (yttria-stabilized  zirconia)  elec¬ 
trodes  with  different  composition,  especially  in  range  from  40:60  to  60:40  wt.%  LSM:YSZ.  The  model  can 
explicitly  capture  the  effects  of  particle  size,  distribution,  and  electrode  composition  on  several  basic 
microstructure  features  and  electrochemical  properties  of  composite  electrodes,  such  as  coordination 
numbers  and  percolation  probability,  total  and  active  three-phase  boundary  length,  and  interfacial  polar¬ 
ization  resistance.  The  model  is  further  used  to  estimate  LSM-YSZ  electrode  performance  with  the  particle 
size  and  distribution  of  the  source  materials.  The  estimation  generally  coincides  with  the  experiment, 
showing  great  potential  in  predicting  power  density  based  on  the  particle  parameters  of  source  materials 
for  SOFCs. 

©  2010  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Solid  oxide  fuel  cell  (SOFC)  is  a  potentially  important  tech¬ 
nique  of  power  sources  due  to  its  high  efficiency  [1].  Nonetheless, 
substantial  technical  questions  remain  challenge,  such  as  how  to 
minimize  the  electrochemical  polarization  resistance  loses  at  their 
electrodes,  especially  the  cathodic  lose  since  activating  oxygen 
molecular  is  more  difficulty  than  hydrogen  [2].  A  typical  cath¬ 
ode  is  a  porous  composite  consisting  of  an  electro-catalyst  phase, 
which  is  also  an  electronic  conductor,  and  an  ionic  conduction 
phase  [3].  Oxygen  reduction  at  the  cathode  is  a  complex  process 
comprising  several  elementary  reactions  and  transport  steps  that 
take  place  at  different  sites  of  the  porous  composite.  Transports 
of  gas-phase  oxygen,  electron,  and  oxygen  ion  are  conducted  at 
the  electrode  pores,  electro-catalyst,  and  ion-conduction  phase, 
respectively.  Adsorbing  oxygen  molecular,  splitting  the  molecular 
to  oxygen  species,  and  diffusing  to  the  charge-transfer  site  occur  at 
the  surface  of  the  electro-catalyst.  And  finally,  the  charge-transfer 
step  is  considered  to  occur  at  sites  where  the  reaction  constituents 
(electrons,  ions,  and  oxygen  species)  can  simultaneously  coexist  [3  ]. 
Lack  of  any  phase,  insufficient  electro-catalyst  surface  area  or  inad¬ 
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equate  site  for  charge-transfer  will  ultimately  lessen  the  overall 
oxygen  reduction  rate,  i.e.  increase  the  polarization  lose. 

Therefore,  the  electrode  properties  such  as  interfacial  polar¬ 
ization  resistance  are  microstructure  dependence.  The  electrode 
microstructure  is  frequently  illustrated  as  a  random  packing  sys¬ 
tem  of  sphere  particles  as  shown  in  Fig.  la.  The  key  microstructure 
parameter  affecting  the  electrochemical  response  is  the  sites  for 
charge-transfer  reaction,  often  termed  as  three-phase  boundaries 
(TPB),  where  ionic  particle,  electronic  particle,  and  pore  contact 
(Fig.  lb).  A  volume  of  work  shows  that  for  many  SOFC  electrodes, 
overall  performance  scales  with  TPB  length  [4,5].  It  should  be  noted 
that  an  active  TPB  is  not  only  a  site  where  the  three  phases  meet 
as  shown  in  Fig.  lb,  but  also  a  site  to  which  the  three  distinct 
reactants  (electrons,  oxygen  ions,  and  oxygen  species)  are  trans¬ 
ported  and  coexist.  While  almost  all  of  pores  are  available  for  gas 
transportation  [6],  not  every  ionic/electronic  particle  is  available 
for  ion/electron  conduction  due  to  stochastic-packing  nature  of 
particles.  That  is,  only  part  of  the  three-phase  meeting  sites  are 
electrochemical  active.  Those  TPBs  are  active  when  their  ionic  and 
electronic  particles  belong  to  percolating  ionic  and  electronic  clus¬ 
ters,  respectively.  Fig.  1  c  shows  a  percolated  cluster.  The  other  TPBs 
such  as  the  insulated  clusters  shown  in  Fig.  Id,  which  are  formed 
with  particles  not  connected  to  the  percolated  ionic/electronic 
phases,  are  not  active  since  either  electrons  or  oxygen  ions  are  not 
available  at  these  sites.  Consequently,  active  TPB  length  is  smaller 
than  the  total  value  and  is  also  structure  sensitive. 
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Nomenclature 

/  probability  density  ( fxnrr 1 ) 

L  length  per  unit  volume  ( pan-2 ) 

LogN  lognormal  distribution 

M  sampling  object  which  can  be  particle  number  or 

volume 

N  normal  distribution 

P  percolation  probability 

R  resistance  ( £2  cm2 ) 

r  particle  radius  (pm) 

S  surface-area  fraction 

U  uniform  distribution 

Z  coordination  number 

Greek  letters 

0  volume  fraction 

a  standard  deviation  of  particle  distribution  (pan) 

fi  mean  particle  radius  (pun) 

0  contact  angle 

p  resistivity  (£2  cm) 

8  electrode  thickness  (pm) 

A  effective  thickness  of  electrode  (pan) 

Subscripts 

i  ionic  particles 

j  electronic  particles 

g  pores 

p  polarization 

TPB  three-phase  boundary 

particles  have  a  distribution 


Theoretically,  the  active  TPB  length  can  be  quantitatively  related 
to  the  two-phase  particle  sizes  and  porosity  [7-9].  However,  ear¬ 
lier  models  assume  that  particles  are  mono-sized,  rather  than  the 
case  with  a  size  distribution  in  question.  Thus,  those  simulations 
do  not  capture  the  effects  of  particle  distributions  on  electrode 
performance.  Recently,  Gokhale  et  al.  [10]  proposed  a  stochastic 
geometry  based  model  (SGM),  relating  total  (active  +  inactive)  TPB 
length  to  mean  size  and  standard  deviation  of  particles.  They  have 
shown  that  mono-sized  particles  lead  to  high  total  TPB  length  com¬ 
pared  to  the  particles  having  size  distributions.  However,  particle 
size  distributions  have  more  complex  effects  on  TPB  length  accord¬ 
ing  to  the  numerical  reconstruction  of  composite  electrodes  [6].  In 
addition,  particle  size  distribution  itself  is  much  more  complicated 
than  a  uniform  distribution  with  standard  deviation.  Thus,  it  is  of 
interest  to  develop  accurate  relationships  between  particle  size  dis¬ 
tributions  and  electrode  properties  such  as  active  TPB  length  and 
interfacial  polarization  resistance. 

In  this  work,  we  propose  an  approach  based  on  percolation  the¬ 
ory  and  particle-layer  model  [9]  to  explicitly  capture  the  effects 
of  particle  size  distributions  on  total  and  active  TPB  lengths  and 
interfacial  polarization  resistance.  Although  the  method  and  results 
presented  here  are  general  and  can  be  applied  to  both  composite 
anodes  and  cathodes,  for  this  work,  we  have  denoted  the  electro¬ 
catalyst  phase,  which  is  also  an  electronic  conduction  phase,  as 
lanthanum  strontium  manganite  (LSM)  and  the  oxygen-ion  con¬ 
ducting  phase  as  yttria-stabilized  zirconia  (YSZ)  to  represent  a 
composite  cathode  and  further  to  validate  our  model. 

The  present  model  shows  that  TPB  length  and  electrochem¬ 
ical  performance  are  very  sensitive  to  particle  sizes  and  their 
distributions.  However,  these  parameters  are  seldom  practically 
known  due  to  the  difficulty  in  phase  identification.  A  combina¬ 
tion  of  focused  ion  beam  (FIB)  [11,12],  Auger  electron  spectroscopy 


Fig.  1.  Schematic  diagrams  of  a  random  parking  system,  (a)  A  porous  two-phase 
composite  electrode  with  particle  size  distributions,  (b)  three-phase  boundary,  (c) 
the  network  of  ionic  particles,  and  (d)  the  electronic  particles  showing  an  isolated 
cluster  that  does  not  contribute  to  active  TPB. 

elemental  mapping,  and  secondary  electron  imaging  has  been  suc¬ 
cessfully  used  for  imaging  the  YSZ  and  LSM  phases  [13].  Recently, 
another  combination  of  chemical  etching  technique  and  atomic 
force  microscopy  is  developed  to  observe  these  phases  [14].  Nev¬ 
ertheless,  revealing  the  size  and  its  distribution  is  very  time 
consuming  and  relays  on  equipments  that  are  not  readily  avail¬ 
able  in  most  research  laboratories.  Here,  an  estimation  method  is 
proposed  to  determine  the  particle  size  and  distribution  from  the 
source  materials. 

The  LSM-YSZ  porous  cathode,  which  is  used  as  the  example 
electrode  in  this  work,  is  usually  processed  with  LSM  and  YSZ 
powders,  typically  by  coating  a  mixture  consisting  of  the  LSM 
and  YSZ  particles  onto  a  dense  YSZ  electrolyte  substrate  using  a 
slurry-based  technique  such  as  screen-printing,  followed  by  co¬ 
sintering  the  two-phase  particles.  The  LSM  and  YSZ  particles  are 
fabricated  with  various  techniques  such  as  solid-state  reaction,  co¬ 
precipitation,  and  combustion.  Whatever,  the  last  step  is  via  firing 
to  form  perovskite  structured  LSM  and  fluorite  structured  YSZ.  The 
firing  temperature  is  usually  close  to  the  co-sintering  temperature, 
which  is  used  to  form  the  composite  cathode.  So,  particle  growth  in 
the  co-sintering  process  is  negligible  when  the  firing  temperature 
is  close  to  or  higher  than  the  co-sintering  temperature.  In  addition, 
in  a  porous  composite  system  consisting  of  two  phases  which  do 
not  react  with  each  other,  particle  growth  (grain  growth)  of  one 
phase  is  limited  by  the  other  phase.  For  example,  compared  with 
pure  LSM  phase,  LSM  particle  growth  in  a  LSM-YSZ  composite  is 
much  less  significant  even  after  it  is  sintered  at  1200°C  [15].  Con¬ 
sequently,  in  most  cases  when  the  sintering  temperature  is  not  too 
high,  the  particle  size  in  the  cathode  is  equal  or  closes  to  that  of 
the  initial  powder.  Thus,  particle  size  and  distribution  of  the  initial 
powder  roughly  represent  the  size  and  distribution  of  the  particle 
in  a  sintered  electrode. 

Accordingly,  the  electrochemical  performance  could  be  pre¬ 
dicted  with  our  model  regarding  the  particle  size  and  its 
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distribution  of  the  initial  LSM  and  YSZ  powders.  So,  in  this  work, 
in  addition  to  general  model  development,  interfacial  polarization 
resistances  of  LSM-YSZ  cathodes,  which  are  prepared  with  differ¬ 
ent  YSZ  particles,  are  theoretically  predicted  and  experimentally 
determined.  The  theoretical  prediction  generally  agrees  with  the 
experimental  results,  suggesting  that  it  is  possible  to  predict  the 
electrode  performance  with  the  original  powder  properties  such 
as  particle  size  and  distribution.  Since  the  electrolyte  conductivity 
is  well  known,  this  work  points  out  a  path  through  which  power 
density  of  SOFCs  (at  least  single  cells)  could  be  predicted  from  the 
characteristics  of  raw  powders.  This  is  a  new  approach  which  is 
here  named  as  a  from-powder-to-power  path. 

2.  Theoretical  development 

Actual  structures  of  composite  electrodes  are  different  from  the 
random  packing  systems,  but  this  approach  has  been  proven  to 
be  successful  in  correlating  micro-structure  parameters  and  elec¬ 
trochemical  performance  of  composite  electrodes  [2,8].  Properties 
of  random  packing  system  of  binary  spherical  particles  can  be 
estimated  using  percolation  theory.  The  concept  of  coordination 
number  is  central  for  the  implementation  of  percolation  theory.  For 
example,  in  a  binary  mixture  of  i  type  (the  ionic  conduction  phase) 
and  j  type  (the  electronic  conduction  phase)  particles  with  radiuses 
of  X[  and  rJt  respectively,  the  coordination  number,  Z2-j,  between  an 
i  type  particle  and  all  j  type  particles  holds  [16], 


Fig.  2.  Definition  of  particle  size  distributions,  (a)  Uniform,  (b)  normal,  (c)  lognor¬ 
mal,  and  (d)  an  actual  distribution  of  YSZ  particles  [17]. 


where  Z  indicates  the  overall  average  coordination  number  of  the 
two-phase  system  and  the  value  is  typically  6;  Sj  is  the  surface-area 
fraction  of  j  type  particles  relative  to  all  particles, 


:  <hhi 
1  <Pi/ri  +  4>j/rj 


(lb) 


where  0;  is  the  volume  fraction  of  phase  i  relative  to  all  solid  phases. 

To  illustrate  the  effects  of  particle  distributions,  uniform  (If), 
normal  (N),  and  lognormal  (LogN)  distributions  are  considered. 
Fig.  2  shows  the  definition  of  these  distributions  with  the  same 
mean  particle  radius,  /x,  and  standard  deviation,  a.  The  probability 
density,/,  is  defined  as, 


m = 


for  r  e  [jii  -  a,  ^  +  cr] 

2(7 


0,  otherwise 
for  the  uniform  distribution  U(/z,cr)  (Fig.  2a); 


f(r)  = 


1 


aV2^A 


exp 


2a2 


for  r  >  0 


(2a) 


(2b) 


for  the  normal  distribution  N(/x,  o )  (Fig.  2b),  where  A  is  normaliza¬ 
tion  constant  making  up  for  the  restriction  of  r  >  0;  and 


m  = 


1 


rnvS 


exp 


-(In  r  -  my 
2  n2 


for  r  >  0 


(2c) 


for  the  lognormal  distribution  Log N(/x,cr)  (Fig.  2c),  where  m  = 
\n(ii2 / y] a2  +  ji2)  and  n  =  yj ln(<r2//x 2  +  1). 

Flowever,  any  type  of  distribution  can  be  implemented,  such 
as  a  real  case  shown  in  Fig.  2d  [17].  In  practical,  SOFC  powders 
always  permit  the  successful  application  of  lognormal  distribution 
[17-19].  As  shown  in  Fig.  2a,  the  differential  element ,/(r)dr,  indi¬ 
cates  the  fraction  of  particles  with  radius  between  r  and  r  +  dr.  It 
should  be  noticed  that  the  physical  meaning  of  the  fraction  can 
be  volume  fraction  or  number  fraction,  represented  as  the  sub¬ 
script  M  =  V  or  N,  respectively.  This  depends  on  sampling  object. 


For  instance,  the  sampling  object  of  a  laser  particle  size  analyzer  is 
particle  volume  [17],  while  the  method  in  Ref.  [6]  and  this  study 
uses  particle  number  as  the  sampling  object. 

Considering  j  type  particles  are  poly-dispersed  with  i  type  par¬ 
ticles  mono-sized,  the  surface-area  fraction  of  j  type  particles  with 
radius  between  y  and  r,  +  dy  relative  to  all  particles  is, 


<Pi¥[ih 

4>i/ri  +  4>j  (l/rj) 


(3a) 


where  <l/ij)  indicates  mean  value  of  1/rj,  defined  as  (\/rf  = 
f  fj/rjdrj.  The  probability  density  function,  /,  is  with  respect  to 
volume  fraction  of  particles.  That  is  to  say  the  sampling  object 
is  particle  volume.  In  principle,  the  probability  density  function 
regarding  particle  number  can  be  derived  from  /as  Bf/r 3,  where 
parameter  B  denotes  normalization  constant.  Clearly,  the  coordi¬ 
nation  number  between  an  i  type  particle  and  j  type  particles  with 
radius  between  r,  and  r,  +  dr,  can  be  estimated  as 


(3b) 


where  underline  represents  particles  have  a  distribution. 

Combining  with  Eqs.  (3a)  and  (3b),  the  coordination  number 
between  an  i  type  particle  and  all  j  type  particles  is, 


zu  = 


dZij  — 


Z  <Pi 
2  4>i/ri  +  0/  (l/rj) 


(4a) 


Considering  i  type  particles  also  have  a  size  distribution /f,  the 
average  coordination  number  between  an  i  type  particle  and  all  j 
type  particles  is  given  by, 


(4b) 
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where  dN2  indicates  the  particle  number  fraction  of  i  type  particles 
with  radius  between  r2  and  r2  +  dr{  relative  to  total  i  type  particle 
numbers, 


dNi  = 


(4c) 


Combining  with  Eqs.  (4a)-(4c),  the  average  coordination  number 
is  obtained  for  the  two-phase  particles  with  size  distributions, 


z.  =?_ i _ 

■-  2  ( W)  <l>i  <i/n)  +  4>j  (Vo) 


fifjdridrj 

(5) 


The  coordination  number  between  i  type  particles,  Z22,  and  j  type 
particles,  Zjj,  can  be  obtained  by  simply  substituting  the  subscripts 
j  and  i  in  Eq.  (5),  respectively. 

Similar  to  binary  system  with  mono-sized  particles  [  7  ] ,  the  prob¬ 
ability  of  an  i  type  particle  belonging  to  the  percolating  cluster  is 
estimated  as, 


Pi  = 


f  4.236 -ZUi 
^  2.472  "" 


(6) 


Consequently,  active  TPB  length  per  unit  volume  can  be  expressed 
as, 


Ltpb  =  PiPj 


2tt  sin  0  min(r2,  rj  )Z2  j 


(JjzfeW i 

4mf/3 


dr, 


(7) 


In  Eq.  (7),  the  integral  item  indicates  total  TPB  length  per  unit  vol¬ 
ume;  ((1  -  0g)0^/47rr3 /3)dr2  is  the  number  of  i  type  particles  with 
radius  between  r2  and  r2  +  dr2  per  unit  volume,  where  (ps  is  electrode 
porosity;  0  indicates  the  contact  angle  between  a  smaller  particle 
and  a  bigger  particle  as  illustrated  in  Fig.  1  b.  The  value  of  0  is  usually 
assumed  between  15°  and  90°  [8,20],  which  is  sensitive  to  experi¬ 
mental  condition.  Combining  Eqs.  (4a)  and  (7)  gives  the  following 
expression  for  active  TPB  length  per  unit  volume: 


brPB  =  ^  sin  OZPiPj 


(1  -  0g)020/ 


0i  ( Vri)  +  0j 


Mj  ff.(  , 

wiC)H  t  h  ,> 


’‘{Tk+w)mdri  m 

It  should  be  noticed  that  when  standard  deviation  of  particles  tends 
to  be  zero,  that  is  the  particles  tend  to  be  mono-sized,  Eq.  (8)  can 
degrade  into  the  formula  based  on  the  classical  percolation  theory 

[91- 

Combining  with  the  particle-layer  model,  interfacial  polariza¬ 
tion  resistance  can  be  estimated  as  [9], 


Rp  =  pfX  coth 


(9a) 


where  8  indicates  electrode  thickness ;  p2  is  resistivity  of  ionic  phase, 

P,  =  P°[(1  -<Pg)4>iPirr  (9b) 

in  which  p9  is  the  intrinsic  resistivity,  and  y  is  the  Bruggeman  factor 
that  accounts  for  tortuous  ion-conduction  pathways  and  the  value 
is  typically  1.5;  A,  is  the  effective  electrode  thickness  within  which 
electrochemical  reactions  mainly  take  place  and  is  estimated  as, 


where  ptpb  is  the  resistivity  of  TPB  derived  from  electrochemical 
reactions. 


3.  Experimental  procedure 

LSM-YSZ  composite  cathodes  were  prepared  onto  YSZ  elec¬ 
trolytes  using  screen-printing,  a  typical  electrode  processing 
technique.  To  prepare  the  YSZ  electrolyte  substrates,  8  mol%  Y203 
stabilized  zirconia  (YSZ)  powder  (Jiujiang  Fanmeiya  Advanced 
Materials  Co.,  Ltd.,  China)  was  cold-pressed  under  20  MPa  to  form 
green  pellets.  The  samples  were  then  sintered  at  1350  °C  for  5h, 
resulting  in  dense  YSZ  substrates  with  the  diameter  of  about  1 0  mm 
and  the  thickness  of  0.5  mm.  LSM  powder  was  synthesized  using  a 
glycine-nitrate  method  [21]  and  fired  at  900  °C  to  form  the  per- 
ovskite  structure.  Commercial  YSZ  powders  with  different  particle 
sizes  were  cooperated  with  LSM  to  form  the  composite  cathodes. 
The  relative  coarse  YSZ  powder  with  big  particle  size  is  represented 
with  C-YSZ,  the  fine  one  with  small  size  is  called  as  F-YSZ,  and  the 
medium  one  is  named  as  M-YSZ.  To  determine  the  particle  size 
and  its  distributions,  the  LSM  and  YSZ  powders  were  ultrasonically 
dispersed  for  5  min  and  observed  with  scanning  electron  micro¬ 
scope  (SEM,  JSM-6700F,  JEOL).  The  particle  size  was  measured  on 
the  SEM  images  using  Matlab®  software.  For  each  sample,  300-400 
particles  were  measured  to  ensure  statistical  accuracy.  The  results 
are  shown  in  Fig.  3.  These  sizes  can  be  well  fitted  to  lognormal 
distribution.  The  mean  particle  radius  and  standard  deviation  are 
also  shown.  To  fabricate  composite  cathodes,  LSM  powder  was 
mixed  with  F-YSZ,  M-YSZ  and  C-YSZ  powders,  respectively,  and 
then  milled  with  ethyl  cellulose  to  form  uniform  slurries.  Symmet¬ 
rical  half-cells  were  prepared  by  coating  the  slurry  on  both  surfaces 
of  the  YSZ  electrolytes  using  the  screen-printing  technique.  The 
cathodes  were  then  co-fired  with  the  electrolytes  at  1000°C  for 
2h  to  produce  a  ~50  p,m-thick  electrode  layer.  Two  parallel- 
samples  were  fabricated  for  each  combination.  Pt  paste  and  Ag 
wires  were  used  for  current  collection.  An  electrochemical  station 
(IM6e,  ZAHNER)  was  used  to  evaluate  the  interfacial  polarization 
resistance.  The  measurement  was  conducted  under  open-circuit 
conditions  with  an  amplitude  10  mV  in  the  frequency  range 
0.01  Hz-1  MHz. 

4.  Results  and  discussion 

4.1.  Coordination  number 

Coordination  number  is  the  core  parameter  of  the  percolation 
theory.  It  is  therefore  discussed  in  this  section  as  a  function  of  par¬ 
ticle  size  distributions  including  distribution  types  and  sampling 
objects.  To  clearly  present  the  distribution  effect  of  one  phase, 
the  other  phase  is  assumed  to  be  monosized.  Let’s  start  with  a 
case  that  j  type  particles  are  not  monosized  but  satisfy  uniform, 
normal  and  lognormal  distributions  with  mean  particle  radius  of 
0.5  p,m  and  standard  deviation  of  Oj.  Their  sampling  objects  can 
be  particle  volume  and  particle  number.  The  electrode  has  a  typi¬ 
cal  composition  of  02  =  4>j  =  0.5.  And  the  i  type  particle  radius  has  a 
monosize  of  0.5  p,m.  Fig.  4  shows  the  effects  of  the  standard  devia¬ 
tion,  Oj  on  coordination  numbers  between  ionic-ionic,  Z22  (Fig.  4a), 
electronic-ionic,  Z2J-  (Fig.  4b),  and  electronic-electronic  particles, 
Zjj  (Fig.  4c).  When  Oj  =  0,  that  is  j  type  particles  are  also  mono-sized 
and  the  radius  is  equal  to  the  other’s,  the  calculated  coordination 
numbers  are  all  3  (Fig.  4).  This  is  in  great  accordance  with  ear¬ 
lier  models  for  composites  consisting  of  mono-dispersed  particles 
[9,16]. 

When  the  sampling  object  is  given,  the  trends  of  coordination 
numbers  are  the  same  for  the  uniform,  normal,  and  lognormal  dis¬ 
tributions.  For  instance,  when  o  increases,  which  means  a  broader 
distribution,  and  the  object  is  particle  number,  Zjj  and  Z2J  decrease 
while  Z2  2  increases.  Increasing  Oj  for  the  number  objected  distribu¬ 
tion  means  the  distribution  is  broadened  almost  equally  (exactly 
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Fig.  3.  SEM  graphs  showing  the  particles  and  distributions  for  LSM,  F-YSZ,  M-YSZ,  and  C-YSZ  powders. 


equal  for  uniform  and  normal  distribution)  by  the  number  of  small 
and  big  particles.  However,  the  big  particles  occupy  more  space. 
Thus,  increasing  Oj  is  equivalent  to  increasing  j  type  particle  size. 
That  is  to  say,  the  larger  particles  act  a  more  important  role  in  the 
whole  particle  range  when  particle  number  is  of  concern.  On  the 
contrary,  the  numbers  exhibit  opposite  trend  when  the  object  is 
volume.  In  this  case,  increasing  a  causes  extraordinary  enlarging 
in  small  particle  number.  So,  decreasing  the  coordination  number 
between  i  particles,  Z22 ,  in  accompanied  by  increasing  the  coordi¬ 
nation  number  related  toj  particles,  i.e.  Zjj  and  Z2J-.  Increasing  in  Z22 
and  Zj  j  indicates  high  percolation  probability  for  i  and  j  phase,  con¬ 
sequently,  the  enhancement  of  ion  and  electron  conducting  ability, 
respectively.  Unfortunately,  as  shown  in  Fig.  4a  and  c,  increasing 
coordination  number  of  one  phase  results  in  the  decrease  of  the 
other.  Practically,  Z22  should  have  a  high  value  since  the  oxygen- 
ion  conductivity  is  much  lower  than  the  electronic  one.  High  Z2J- 
value  indicates  high  TPB  length,  therefore  improving  electrochem¬ 
ical  performance. 

Fig.  4  also  shows  that,  although  the  curvilinear  trend  is  the  same 
for  the  distributions  termed  with  uniform,  normal,  and  lognor¬ 
mal,  absolute  change  in  coordination  number  is  quite  difference. 
While  the  changes  are  comparable  for  the  normal  and  lognormal 
distributions,  the  change  for  the  uniform  distribution  is  much  less 
significant  even  though  a  is  the  same.  This  is  because  the  ranges  of 
the  normal  and  lognormal  distributions  are  comparable  and  they 
are  much  wider  than  that  of  the  uniform  distribution,  as  shown  in 
Fig.  2. 


4.2.  Percolation  probability 

Fig.  5  shows  the  percolation  probability,  Pj,  as  a  function  of  vol¬ 
ume  fraction,  </>j,  when  j  particles  have  distributions  with  mean 
radius  of  0.5  |xm  and  standard  deviation  of  0.2  |xm  and  i  particles  are 
monosized,  r2  =  0.5  [xm.  When  both  phases  have  the  same  size,  i.e. 
r,  =  0.5  [xm  as  shown  in  Fig.  5,  the  same  curve  can  be  also  calculated 
using  our  previous  model  for  composites  composed  of  monosized 
particles  [9].  Compared  with  the  mono-sized  case,  when  the  sam¬ 
pling  object  is  particle  number,  percolation  probability  decreases 
and  volume  fraction  threshold  shifts  to  a  high  level.  The  predic¬ 
tion  suggests  more  j  type  phase  should  be  used  to  maintain  the 
same  percolation  probability.  This  is  in  accordance  with  the  results 
derived  from  3D  numerical  simulation  [6].  On  the  other  hand, 
when  the  sampling  object  is  particle  volume,  percolation  proba¬ 
bility  increases,  and  volume  fraction  threshold  tends  to  be  low. 
The  tendency  is  much  more  significant  for  the  normal  and  lognor¬ 
mal  distributions  because  of  the  wider  range  compared  with  the 
uniform  distribution. 


4.3.  Active  TPB  length 

Fig.  6  shows  how  the  particle  size  distribution  can  affect  active 
TPB  length  at  various  electrode  compositions.  Normal  distribution 
is  implemented  to  compare  with  the  3D  numerical  reconstruction 
results  [6]  and  SGM  approach  [  1 0].  The  contact  angle  6  is  assumed  as 
34.4°  according  to  the  cross-section  of  3D  reconstructed  electrode. 
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Fig.  4.  Effects  of  standard  deviation,  Oj,  on  coordination  numbers  for/  type  particles  with  uniform  ( U ),  normal  (N),  and  lognormal  (LogN)  distributions  (0,-  =  0j  =  0.5,  r,-  =  0.5  (Jim, 
lij  =  0.5  |jim).  (a)  Coordination  number  within  i  type  particles,  (b)  between  i  type  and  j  type  particles,  and  (c)  within  j  type  particles. 


Fig.  6a  shows  the  comparison  of  active  TPB  length  derived  from 
3D  reconstruction,  SGM  model,  and  the  present  approach.  When 
the  sampling  object  is  number,  all  the  three  results  show  that  the 
active  TPB  length  derived  with  particle  distributions  of  ri  =  rj  =  NN 
(0.25  |jim,  0.1  [Jim)  is  lower  than  that  with  monosized  particles, 
r2  =  rj  =  0.25  pm.  This  phenomenon  is  mainly  caused  by  the  decrease 
of  coordination  number  Z2J-  (Fig.  4b)  and  percolation  probability 
(Fig.  5)  when  the  particles  have  a  distribution  termed  with  num¬ 
ber.  These  methods  also  demonstrate  that  the  highest  TPB  length 
is  observed  at  </>2  =  </>j  when  the  two  phases  have  the  same  particle 


Fig.  5.  Effects  of  /  phase  volume  fraction,  0/,  on  its  percolation  probability  for  various 
j  type  particle  distributions  with  /ij  =  0.5  p,m  and  (Tj  =  0.2  |jim.  The  i  type  phase  is 
monosized,  r,  =  0.5  p,m. 


sizes  and  distributions.  Flowever,  there  are  differences  in  the  results 
derived  by  the  three  methods.  Our  model  calculation  agrees  well 
with  the  3D  reconstruction  result  in  the  whole  composition  range 
limited  by  the  percolation  requirement.  Compared  with  the  3D 
reconstruction  method,  the  model  can  be  easily  conducted  because 
it  is  analytic.  Indeed,  by  given  specific  knowledge  for  particular 
electrode  micro-structures,  the  model  can  be  applied  to  predict 
general  results  and  explicitly  capture  the  dependence  of  active 
TPB  length  on  particle  size  distribution.  Fig.  6a  also  shows  that  our 
model  produces  better  agreement  with  3D  reconstruction  data  than 
SGM  model  in  a  wider  composition  range.  To  the  best  of  our  knowl¬ 
edge,  SGM  approach  is  the  only  available  model  considering  particle 
size  distributions  on  TPB  length  calculation.  Nonetheless,  SGM  cal¬ 
culates  total  TPB  length  including  effective  and  non-effective  sites. 
As  shown  in  Fig.  1 ,  not  all  the  TPBs  are  active  for  electrochemical 
reaction  since  some  sites  are  only  physical  interfaces  where  the  ion- 
conduction  phase,  electro-catalyst  phase,  and  pore  meet,  not  those 
sites  for  charge-transfer  reaction  where  oxygen  ions,  electrons,  and 
oxygen  molecular  coexist.  It  should  be  noted  that  at  typical  com¬ 
position  range  0j  =  O.4-O.6,  about  90%  of  total  TPBs  is  active  as 
shown  in  Fig.  5.  Furthermore,  SGM  concludes  that  the  maximum 
TPB  length  always  corresponds  to  the  composition  of  02  =  =  0.5 

disregarding  difference  of  particle  size  distribution.  It  seems  cor¬ 
rect  as  suggested  by  Fig.  6a.  Flowever,  the  3D  reconstruction  and 
our  model  show  that  the  composition  corresponding  to  the  maxi¬ 
mum  TPB  length  is  a  function  of  particle  size  and  its  distribution, 
which  is  discussed  in  the  following  sections. 

It  is  of  interest  to  note  that  if  the  particle  radiuses  of  both 
ionic  and  electronic  phases  are  the  same,  and  they  have  the  same 
distributions,  the  volume  fraction  of  ionic  or  electronic  phase 
corresponding  to  the  maximum  active  TPB  length  will  be  0.5.  How- 
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Fig.  6.  TPB  length  for  composite  electrodes  consisting  of  various  particles.  Calculation  is  conducted  based  on  0g  =  0.3  and  6  =  34.4°  according  to  3D  reconstruction  results,  (a) 
Total  and  active;  (b)  active  TPB  length  derived  by  3D  reconstruction  [6],  SGM  model  (total  TPB  length)  [10],  and  the  present  approach  for  composite  electrodes  composed 
of  monosized  particles  (rj  =  r,  =  0.25  p,m)  and  of  particles  with  normal  size  distribution,  Nn  (0.25  p,m,  0.1  |xm).  The  sampling  object  is  particle  number;  (c)  active  TPB  length 
derived  from  different  distributions  sampled  with  particle  volume;  and  (d)  effects  of  standard  deviation  and  sampling  objects  on  active  TPB  length. 


ever,  the  volume  fraction  can  be  different  when  the  mean  particle 
radiuses  of  both  phases  are  same  and  particle  size  distributions  are 
different.  For  instance,  in  case  of  r2  =  0.25  [xm  and  rj  =  NN  (0.25  pun, 
0.1  |xm),  Fig.  6b  demonstrates  that  (f>j  corresponding  to  the  maxi¬ 
mum  active  TPB  length  shifts  to  0.55  from  0.5  for  r2  =  r,  =  0.25  |xm. 
This  is  caused  by  the  increase  of  (j)j  threshold  as  shown  in  Fig.  5  and 
the  decrease  of  02  threshold. 

When  sampling  method  is  particle  volume,  the  results  are  dif¬ 
ferent.  For  instance,  as  shown  in  Fig.  6c,  in  the  case  of  r2  =  0.25  p,m 
and  rj  =  Ny  (0.25  |xm,  0.1  |xm),  the  volume  fraction,  <pj,  relating  to  the 
maximum  TPB  length  shifts  to  0.44,  being  contrary  to  that  shown  in 
Fig.  6b.  Reasons  are  the  increase  of  Pj  (Fig.  5)  and  the  decrease  of  P2 
when  the  sampling  method  changes  from  particle  number  to  vol¬ 
ume.  On  the  other  hand,  in  the  case  of  r2  =  r,  =  Nv  (0.25  |xm,  0.1  pan), 
Fig.  6c  shows  that  the  active  TPB  length  is  remarkably  enhanced 
compared  with  the  mono-sized  case,  which  is  also  different  with 
the  results  shown  in  Fig.  6a.  This  is  the  result  ofZ2J-  increase  (Fig.  4b). 
Combining  Fig.  6c  with  Fig.  6a,  it  can  be  concluded  that  when  i  and 
j  type  particles  have  the  same  distribution  with  respect  to  particle 
volume/number,  the  active  TPB  length  can  increase/decrease  com¬ 
paring  to  its  mean  size  case.  This  is  further  confirmed  in  Fig.  6d. 
However,  it  should  be  noted  that  when  o  is  small,  there  is  no  obvi¬ 
ous  difference  between  the  two  types  of  particle  size  distributions, 
and  smaller  number  of  total  particles  is  required  to  fill  the  electrode 
domain  by  comparing  with  large  poly-dispersed  case.  Therefore, 
active  TPB  length  is  reduced. 


4.4.  Application  to  a  real  electrode  system 

In  a  real  electrode  system,  particle  size  and  distribution  are 
rarely  the  same  for  the  two  phases.  Thus  it  is  of  great  interest  to 
apply  our  approach  to  a  real  system  rather  than  the  3D  numeri¬ 
cal  reconstruction.  But  data  is  seldom  reported  on  the  structural 
details  of  composite  electrodes.  Recently,  Barnett  et  al.  report  the 


microstructural  features  of  LSM-YSZ  electrodes  with  compositions 
ranging  from  30:70  to  70:30 wt.%  LSM:YSZ,  including  total  and 
active  TPB  length,  which  are  determined  using  FIB-SEM  technique 
[12].  The  electrodes  are  prepared  with  the  screen-printing  method 
and  sintered  at  1175  °C  for  1  h.  Before  printing,  LSM  and  YSZ  pow¬ 
ders  are  mixed  through  ball-milling.  Data  is  not  directly  available 
on  the  particle  size  distribution  for  the  LSM  and  YSZ  particles  in 
the  electrodes.  However,  using  the  intercept  method  [22],  it  can 
be  estimated  from  the  dispersed  particles  on  the  segmented  2D 
images,  which  are  shown  in  Fig.  3  of  Ref.  [12].  Fig.  7a  shows  the 
results  for  particle  radius  and  its  distribution.  Clearly,  LSM  and  YSZ 
particles  are  not  of  the  same  size.  They  can  be  represented  with  log¬ 
normal  distributions,  Log  %  (0.27  |xm,  0.08  |xm)  for  LSM  and  Log  Nn 
(0.1 9  |xm,  0.05  |xm)  for  YSZ.  The  total  and  active  TPB  length  could  be 
calculated  using  the  present  approach.  According  to  the  suggestion 
of  Ref.  [23],  contact  angle  0  is  assumed  as  19.5°  for  the  LSM-YSZ 
systems.  Fig.  7b  shows  the  comparison  between  the  prediction  and 
the  experimental  results  obtained  with  FIB-SEM  by  Barnett  et  al. 
The  TPB  length  is  also  calculated  with  SGM  model,  which  deals 
only  with  the  total  length.  Our  calculation  results  agree  well  with 
experimental  data  in  both  the  active  and  the  total  TPB  lengths. 

However,  it  should  be  noticed  that  the  present  model  predicts 
zero  active  TPB  length  beyond  the  composition  threshold.  As  model 
calculation  shown  in  Fig.  7c,  when  YSZ  composition  is  beyond  the 
threshold  (0ysz>O-61),  LSM  particles  cannot  form  a  percolating 
cluster,  represented  by  zero  percolation  probability.  Thus  there  is 
no  active  TPB  theoretically  exists  within  the  electrode,  even  almost 
all  YSZ  particles  are  percolated.  But  there  exists  a  small  amount 
of  active  TPBs  according  to  the  FIB-SEM  results.  These  TPBs  are 
mainly  close  to  the  current  collecting  layer.  They  can  be  formed  by 
the  isolated  LSM  clusters  (Fig.  Id)  contacting  with  the  current  col¬ 
lector  and  the  percolating  YSZ  particles.  The  isolated  clusters  are 
not  theoretically  considered  as  the  active  sites  in  the  percolation 
theory.  In  addition,  the  finite  domain  size  used  for  FIB-SEM  may 
be  too  small,  throughout  which  the  probability  of  LSM  particles 


1990 


Y.  Zhang  et  al.  /  Journal  of  Power  Sources  196(2011)  1983-1991 


Fig.  7.  Comparison  between  model  prediction  and  experimental  results  for  LSM-YSZ  electrodes  [12].  (a)  Particle  size  distribution  of  LSM  and  YSZ  particles  derived  from  Fig. 
3  of  Ref.  [12];  (b)  TPB  length  determined  with  present  approach,  SGM  model  and  FIB-SEM  method;  (c)  percolation  probabilities;  and  (d)  interfacial  polarization  resistance 
(0g  =  O.5). 


can  be  overestimated.  Nevertheless,  this  difference  suggests  that 
those  isolated  particles  and  clusters,  which  are  close  to  the  current- 
collector/cathode  interface  or  electrolyte/cathode  interface,  should 
be  active  for  electron  or  oxygen-ion  conduction  and  thus  have  to 
be  considered  in  the  future  work. 

Fig.  7d  shows  the  comparison  of  interfacial  polarization 
resistances  obtained  with  model  prediction  and  experimental  mea¬ 
surements.  The  resistivities  of  YSZ  and  TPB  are  derived  from  Ref. 
[9].  Near  the  composition  of  0Y sz  =  0-5,  the  prediction  seems  to  be  in 
good  agreement  with  the  experimental  data.  For  the  electrode  with 
porosity  of  0.5  and  thickness  of  11  p,m  [12],  the  predicted  interfa¬ 
cial  resistance  is  0.50  £2  cm2  at  800  °C,  being  comparable  with  the 
reported  0.51  £2  cm2  for  0Ysz  =  O-5.  When  0YSz>O.6  or  0YSz<O.4, 
the  predicted  resistance  is  not  consistence  with  the  experimental 
data.  When  LSM  content  is  low  (0YSz  >  0.6),  the  predicted  resistance 
is  higher  than  the  measured  one.  At  this  composition,  almost  all 
the  YSZ  particles  belong  to  a  huge  percolated  YSZ  cluster  filling 
all  the  electrode  space,  while  almost  all  the  LSM  particles  are  not 
percolated  but  form  some  isolated  clusters.  TPB  formed  with  the 
isolated  clusters  are  possibly  electrochemical  active  when  the  clus¬ 
ters  are  on  the  current  collector  side.  For  these  TPBs,  oxygen  ions 
are  readily  available  through  the  percolated  YSZ  phase  and  elec¬ 
trons  are  also  available  through  the  isolated  clusters,  which  connect 
to  the  current  collector.  This  underestimates  the  electrochemi¬ 
cal  activity  when  LSM  content  is  low.  On  the  other  hand,  when 
YSZ  content  is  lower  (0YSz  <  0-4),  the  predicted  resistance  is  lower 
than  the  measured  data.  In  this  case,  electrons  are  readily  available 
throughout  the  whole  cathode  since  LSM  percolation  probability 
is  close  to  100%.  Fortunately,  YSZ  particles  are  also  percolated  and 
the  percolation  probability  is  very  high  in  the  experimental  com¬ 
position  range  due  to  the  relatively  small  particles  (Fig.  7a  and 
c).  So,  TPB  formed  with  percolated  YSZ  particles  dominates  the 
electrochemical  activity  although  TPB  created  with  isolated  YSZ 
clusters,  which  connect  to  the  electrolyte  layer,  contributes  also 


to  the  activity.  The  resistance  is  calculated  by  combining  the  active 
TPB  length  with  particle-layer  model  [9],  in  which  Bruggeman  fac¬ 
tor  is  assumed  to  be  a  constant.  The  factor  is  a  function  of  the 
tortuosity  of  ion-conduction  paths,  increasing  with  the  decrease 
of  percolation  probability.  Experimental  results  also  show  that  the 
tortuosity  of  percolated  YSZ  increased  rapidly  with  the  decrease 
of  YSZ  content  [12].  The  increased  tortuosity  infers  the  enlarged 
distance  for  oxygen  ion  transportation,  leading  to  the  discrepancy 
between  the  model  prediction  and  the  experimental  data.  Although 
discrepancies  exist  when  the  volume  fractions  are  quite  different 
for  the  two  phases,  the  model  is  useful  in  predicting  electrochem¬ 
ical  performance  of  practical  electrodes  since  in  these  systems  the 
two  phase  contents  are  quite  close  to  achieve  high  electrode  activ¬ 
ities. 

4.5.  TPB  length  and  interfacial  polarization  resistance  estimated 
with  powder  sizes 

Since  FIB-SEM  instruments  are  not  readily  available  for  nor¬ 
mal  laboratories,  it  is  worthy  to  use  the  particle  size  of  the  source 
powders  to  predict  TPB  length  as  well  as  interfacial  polarization 
resistance  (area  specific  resistance),  especially  when  the  electrodes 
are  formed  at  relatively  low  sintering  temperatures.  YSZ  powders  of 
different  size  are  chosen  for  such  evaluation  due  to  the  fact  that  it  is 
often  much  harder  to  sinter  YSZ  than  LSM,  particularly  at  the  same 
sintering  temperature.  Fig.  8  shows  the  estimated  active  TPB  length 
for  F-YSZ-LSM,  M-YSZ-LSM  and  C-YSZ-LSM  electrodes.  The  pow¬ 
der  particle  distributions  are  fitted  using  lognormal  distribution  as 
shown  in  Fig.  2.  The  sampling  object  is  particle  number.  And  the 
contact  angle  is  assumed  as  19.5°.  Generally,  TPB  length  decreases 
when  YSZ  size  increases.  The  YSZ  volume  fraction  corresponding  to 
the  maximal  active  TPB  length  is  small  for  fine  YSZ.  This  is  caused 
mainly  by  the  smaller  particles  with  smaller  standard  deviation  of 
F-YSZ  than  that  of  LSM  (Fig.  2). 
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Fig.  8.  Active  TPB  length  for  F-YSZ-LSM,  M-YSZ-LSM  and  C-YSZ-LSM  electrodes 
as  a  function  ofYSZ  loading  (0g  =  O.3)  calculated  with  particle  size  from  the  source 
powders. 


Fig.  9.  Interfacial  polarization  resistances  of  F-YSZ-LSM,  M-YSZ-LSM  and  C- 
YSZ-LSM  electrodes  as  a  function  of  temperature  (</>g  =  0.3,  <5  =  50|jim).  Data  is 
estimated  with  size  parameters  of  the  source  powders. 

Fig.  9  shows  the  comparison  in  interfacial  polarization  resis¬ 
tances  as  estimated  from  the  powder  size  and  determined  with 
impedance  spectroscopy  for  F-YSZ-LSM,  M-YSZ-LSM,  and  C- 
YSZ-LSM  electrodes.  YSZ  volume  fraction  is  0.26,  0.48  and  0.64 
for  F-YSZ-LSM,  M-YSZ-LSM,  and  C-YSZ-LSM  electrodes,  respec¬ 
tively.  The  composition  is  set  according  to  the  maximum  TPB  length 
(Fig.  8).  The  resistivity  is  a  function  of  temperature,  which  is  derived 
from  Ref.  [9].  As  shown  in  Fig.  9,  the  experimental  results  gener¬ 
ally  agree  with  the  estimation,  especially  for  the  electrode  based  on 
fine  YSZ.  For  example,  the  experimental  resistance  of  F-YSZ-LSM 
electrode  is  0.26  £2  cm2  at  800  °C,  while  the  estimated  resistance 
is  0.25  £2  cm2.  For  F-YSZ-LSM  electrode,  the  average  difference 
between  the  experimental  and  estimated  results  is  only  3%.  FIow- 
ever,  the  difference  is  big  for  the  electrode  prepared  with  the  coarse 
YSZ,  -18%.  But  the  prediction  is  also  not  too  bad.  Estimation  with 
source  particles  is  limited  by  the  sintering  process,  in  which  sin¬ 
tering  temperature,  time,  and  even  heating  rate  affect  the  size  and 
shape  of  the  particles  in  the  resulted  electrodes.  Such  estimation 
will  become  much  more  reliable  and  can  reveal  the  real  electrode 
performance  if  the  sintering  process  is  clear,  especially  the  process 
for  electro-catalyst,  which  is  much  easier  to  be  sintered  than  the 
electrolyte. 


5.  Conclusion 

An  analytical  model  is  proposed  to  estimate  the  effects  of  parti¬ 
cle  size  and  its  distribution  on  the  electrochemical  performance  and 
microstructure  features  of  SOFC  electrodes.  Active  TPB  length  per 
unit  volume  and  interfacial  polarization  resistance  are  calculated 
and  the  results  are  in  accordance  with  literature  and  experimental 
data.  The  model  produces  explicit  insight  and  quantitative  design 
guidance  for  electrode  fabrication,  and  leads  to  the  following  con¬ 
clusions: 

(1)  The  effects  of  particle  size  distribution  are  sensitive  to  its 
sampling  object.  When  the  sampling  object  is  particle  vol¬ 
ume/number,  the  particles  with  a  wide  size  distribution  lead  to 
a  higher/lower  maximum  TPB  length  as  compared  to  its  mono¬ 
size  particles.  The  difference  is  more  significant  for  a  wider 
distribution. 

(2)  Distribution  with  the  sampling  object  as  volume/number  can 
shift  the  optimal  composition  corresponding  to  the  maximum 
active  TPB  length.  When  j  type  particles  have  a  wider  distribu¬ 
tion  compared  with  i  type  particles,  less/more  j  phase  loading 
is  required  compared  with  its  mono-size  case. 

(3)  It  is  possible  to  estimate  electrode  performance  with  the  par¬ 
ticle  size  and  distribution  of  source  materials.  Consequently, 
SOFC  performance  can  be  theoretically  predicted  with  size 
parameters  of  the  starting  oxide  powders,  suggesting  a  from- 
powder-to-power  approach. 
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